Enhanced Deflection of Light Ray by Atomic Ensemble on Coherent Population 

Oscillation 
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In recent experiments [e.g., Nature Physics 2, 332 (2006)], the enhanced light deflection in an 
atomic ensemble due to inhomogeneous fields is demonstrated by the electromagnetically induced 
transparency (EIT) based mechanism. In this paper, we explore an different mechanism for the 
similar phenomenon of the enhanced light deflection. This mechanism is based on the coherent 
population oscillation, which leads to the hole burning in the absorption spectrum. The medium 
causing the deflection of probe light is an ensemble of two-level atoms manipulated by a strong 
controlled field on the two photon resonances. In the large detuning condition, the response of 
the medium to the pump field and signal field is obtained with steady state approximation. And 
it is found that after the probe field travels across the medium, the signal ray bends due to the 
spatial-dependent profile of the control beam. 

PACS numbers: 42.50.Md, 03.65.Sq, 42.25.Bs 



I. INTRODUCTION 

Recently, slow light propagation has attracted a great 
deal of attention due to the fundamental aspects of non- 
trivial light pulse manipulation [l[ and possible applica- 
tions for optical delay lines [3, f° r quantum computing 
and quantum communication H, and for developing sen- 
sitive measurement techniques Q. The experiments for 
slow light are carried out on various types of materials 
such as cold sodium atoms [E H, a tom vapors fj\. M. [qL fiol] . 
and solids state system [lllll2ul3l. Ibl. [Ta| . 

The physical mechanism of slow light is to create a 
narrow transparency window (burn a hole) within the 
absorption line by an intense coupling laser field. Also 
a sharp normal dispersion of the refractive index is ac- 
companied within this narrow window. Thus both a 
low group velocity of light and an enhanced transmis- 
sion of light can be achieved. A basic way to produce 
slow light is electromagnetically induced transparency 
(EIT) [E S H El. The EIT effect dually happens in 
the so-called A-type atomic system, which contains two 
lower states with separate couplings to an excited state 
through two electromagnetic fields (probe and control 
light). When the absorption of light by both transitions is 
suppressed due to destructive interference between exci- 
tation pathways to the upper level, the medium becomes 
transparent with respect to the probe field. The group 
velocity of the light depends on the parameters of the 
control field [U,[20j]. 

Another method to producing slow light is the coherent 
population oscillation (CPO) @, EE E3, EE E3 • In tnis 
sense, the transition is excited by a probe field and a con- 
trol light from the ground state to the excited state of a 
two-level system. The beating due to the slight detuning 
between pump and probe lights leads to a periodic mod- 
ulation of the atomic population, which create a narrow 
spectral dip in the probe absorptive spectrum. Due to its 
insensitive to the dephasing of the atomic coherence, the 



slow light has been found at room temperature [2|], and 
its speed can be reduced as low as a few tens of meters 
per second (l3| . 

Actually, the conventional studies of slow-light phe- 
nomenon focus on various effects in the time-domain, 
but most recently, much attention has been paid to light 
propagation in the spatial domain. The effect of an ex- 
ternal field with spatially inhomogeneous distribution has 
been studied [H, H Et HI . It is found that the light 
ray bends when a magnetic field with small gradient ver- 
tical to the propagation direction is applied to a EIT 
medium jT^. Also in a EIT medium, the light deflec- 
tion is explicitly observed when the atomic ensemble is 
driven by a optical field with an inhomogeneous profile 
along transverse direction [2l| . 

For the slow light phenomena due to CPO, almost all 
theoretical treatments are confined to the homogeneous 
transverse spatial distribution of the control field. In this 
paper, we address how the spatial profile of the control 
field along transverse direction affects the propagation of 
the probe beam for the CPO based slow light . Here we 
treat the external fields classically. After investigating 
the atomic response to the control field and the probe 
field by the perturbation approach, the spatial motions 
of the control field is governed by an effective nonlin- 
ear Schrodinger equation in the low intensity limit. It is 
found that the control field can propagate for long dis- 
tances with an invariant transverse profile. And with 
an effective potential induced by the steady atomic re- 
sponse, the spatial motions of the probe field is described 
by an effective Schrodinger equation in the limits of lin- 
ear response. And therefore the deflection of the light 
ray comes out straightforwardly. 

This paper is organized as follow: in sec. [ill we present 
the theoretical model for a two-level atomic ensemble in- 
teracting with a control field and a much weaker probe 
field. In Sec. IIII1 the perturbation theory is applied to 
obtain the atomic motion equation which is related to 
the response to the external field. In Sec. IIV1 we derive 
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FIG. 1: (a) Atomic two- level system interacting with the 
probe light and a control light. Qi denotes the Rabi frequency 
of the control and probe field, (b) Configuration of the optical 
beams inside the atomic medium. 



a system of equations which govern the spatial motion of 
the control field and the probe field in the optically con- 
trolled atomic medium. Then the deflection of the probe 
beam is investigated by the transverse spatial profile of 
the control light. In Sec.|Vl we make our conclusion. 



II. OPTICALLY CONTROLLED TWO LEVEL 
ATOMS 



The system for producing slow light by CPO involves 
three subsystems: the atomic medium, a control field and 
a weaker signal field. The atomic medium is an ensemble 
of N identical atoms confined in a cell ABCD shown in 
Fig. [IJb). Each of atoms has two internal states - the 
ground state \g) and the excited state |e). The atomic 
transition is driven by two optical fields, a weaker probe 
field and a stronger control field, as shown in the top 
panel of Fig.[TJ The energy difference between \g) and |e) 
is denoted by u> eg — uj e — to g . A = v c — uj eg is the detuning 
between the atomic transition and the control field with 
carried frequency v c . In this literature, large detuning 
is considered. The probe field is slightly detuned from 
the control field, and it is denoted by d = v c — v p . The 
detuning between the control and probe light leads to 
periodic modulation of the two-level system population 
oscillation, which substantially reduces the group veloc- 

As shown in Fig. [IJb), both the probe field and control 
field propagate parallel along the z-direction with wave 
number k and k c respectively. The Hamiltonian 



(1) 



medium, we introduce the collective atomic operator Q 

^M) = ^^<,W, (2) 



N r 

rj£N 

which average over a small but macroscopic volume V 
containing many atoms N r = (N/V) dr»l around po- 
sition r, where N is the total number of atoms. Then the 
unperturbed atomic Hamiltonian reads 



H (A) 



N 
V 



d 3 r(uj e a ee + u g a gg ) . 



(3) 



where we have neglected the kinetic term of atoms. The 
interaction Hamiltonian describes the coupling be- 
tween the external field and the atomic ensemble. Under 
the electric-dipole approximation and the rotating-wave 
approximation, the interaction Hamiltonian H^ 1 ' reads 



H {i) 



N 
V 



d 3 r ( fL,<7. 



p^eg 



h.c. 



where 



Qj = d eg Ej /h, j =p,c 



(4) 



(5) 



are the Rabi frequencies associated with the correspond- 
ing external fields. is the positive frequency part of 
the corresponding external field. 

By introducing the slow varying variables for both 
fields 



E+ = E c (r)e i ^ z -^ t \ 
and for atomic transition operator 



(6a) 
(6b) 

(7) 



the dynamics of this system is described by the interac- 
tion Hamiltonian 



H, 



N 
V 
+h.c 



n 



(E p( 



i(5t~k$ z) 



EAn 



(8) 



is composed of two terms. To describe the atomic 



in the rotating reference frame, where k$ — k c — k p and 
d eg is dipole transition element. Here we have assumed 
that the control beam is monochromatic wave, however 
the slow varying variable E p for the probe light is a su- 
perposition of monochromatic waves, which have a small 
frequency variance Av around the mean frequency u p . 



III. ATOMIC RESPONSE 

It is well-known that, when atoms are subjected to 
an electric field, the applied field displaces the posi- 
tive charges and the negative charges in atoms from 
their usual positions. This small movement that positive 
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charges in one direction and negative ones in the other 
will result in collective induced electric-dipole moments. 
Every dipole behaves collectively to give a response to 
the light. We now investigate the response of the two 
level atomic ensemble driven by a monochromatic clas- 
sical field. The equations for the atomic coherence are 
described by the density matrix elements 



i-T^-VenO-ne + h.C. 



VZ,a ea + h.c. 



h e9 eg 

iAa 



d* 

l -fr V eg {Vee ~ <Jgg) ■ 



where 



V eg = E+e^ ksi 



+ E+ 



(9a) 
(9b) 
(9c) 

(10) 



is the amplitude of the applied optical field. We con- 
sider the ensemble of closed two-level atoms. Due to the 
spontaneous emission, the population of the upper level 
will decay. We assume that the upper level \e) decays to 
the lower level |<?) at a rate 71, which means the lifetime 
of the upper level is given by T\ = 7^ . As the two- 
level is closed, any population that leaves the upper level 
enters the ground level. Thus T\ is the life time of the 
population difference of the ground \g) and excited |e) 
states. We also assume that the atomic dipole moment 
dephases at a rate 72 = T 2 ~ . By adding decay terms phe- 
nomenologically, the dynamic of this atomic ensemble is 
described by 

2d* 

W = -71 (w - W eq ) + 1—^-VegOge + h.C. (11a) 

d* 

a eg = -(iA + "/ 2 )a eg -i-^V eg w. (lib) 
where w is the population inversion operator defined as 



'39' 



(12) 



<(:'</ 



is the population difference in thermal equilibrium. 
Eq. (TTTj) cannot readily be solved exactly. As the probe 
field is much weaker than the control field, we seek a so- 
lution to the equations of the atomic motion that is cor- 
rect to all orders in the amplitude of the control field and 
is correct to lowest order in the amplitude of the probe 
field. Though the polarization induced by the external 
fields in the ensemble of two-level atoms is described by 
the equation 



P (r, t) 



V ^ ^ ^9 



a i{kjZ-Vjt) 



which generally includes components not only at the field 
frequencies v c and v p but also at v p ± k {y c — u p ) with k 
integral, only the frequencies v c , v p and v c + S will occur 
in the atomic response to lowest order of the amplitude 



E p . Thus we treat the atomic equations perturbatively, 
and require that the steady-state solution of Eq. (jTTj) be 
of the form H 113 



w = w w +w M e m*-*t) + w (-) e -^ z - 5t )(i4b) 

where i// ' and (Teg are the population and polarization 
in the absence of the probe field, and other terms are in 
the first order of E p . Actually, E qs. (fT4")) is a truncation 
of the Floquet basis expansion [24]]. The dynamics of 
population and dipole moment of the atomic ensemble 
are given by the equations in the zeroth order of E p 



:,(0) _ 



■ 7 (0) = _ 

eg 
v(o) 



_ 7lW (0) + JlWeg + 2i (p*a$ - h.c}) (15a) 
(zA- 72 )aW +m c wM 



and the first order of E p 

w M = (i<5- 7 i)w (_) ~i2Q p a { e °J 

+i2 (Cltf-) - h.c.) 
&tf = [i(S-A)- l2 ]a^-iW c 



,(-) 



ge 



+iil p w ( > ' > + ifi c io (_) 



where we have introduce the Rabi frequencies 



tt 3 = —Ej, j = c,p. 



(15b) 
(15c) 

(16a) 

(16b) 
(16c) 



(17) 



Under the adiabatic approximation that the evolution 
of the atomic system is much faster than the temporal 
change of the radiation field, we set the left hand sides 
of Eqs. (fTS"! and (fT6|) to zero. Therefore, the population 
and dipole moment in the zeroth order of E p read 



,(0) 



t<°> = 
ge 



7i (A 2 + 7l) Weq 
7i (A 2 + 7 2 2 ) + 47 2 |fi c 
7i («7a - A) O c 



7i (A 2 



•7l) 



-4 72 |ft c 



7 We 



(18a) 
(18b) 



and the dipole moment in the first order of the probe 
field amplitude reads 



D (i 7 2 - A) 



D («7 2 - A) (5 + A + i l2 ) 
2\VL C \ 2 {8 + 2i l2 ){5- A + i l2 ) 

\ 

D (i 72 - A) (6 + A + ij 2 ) 



/ '0,, (19) 



(13) where 



D 



(<5 + i7i)0*- A + i l2 ){8 + A + i l2 ) (20) 

-4\n c \ 2 (i l2 + s) 

Here a g J and a ge •* determine the atomic response to the 
control and probe fields respectively. 
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IV. PROPAGATION OF LIGHTS IN AN 
TWO-LEVEL ATOMIC MEDIUM 

We now consider the deflection of the probe light in the 
rectangular medium due to the transverse spatial profile 
of the control light shown in Fig. For a nonmag- 

netic medium with no free charges and no free currents, 
a driven wave equation 



V^-eoMog^MogrP 



(21) 



can be obtained from the Maxwell equations, where V 2 
is the Laplacian operator, and E is the sum of ampli- 
tudes of the control and probe beam. Here the dielectric 
response P acts as an effective source to generate the 
electromagnetic field. In the paraxial approximation and 
slowing varying amplitude approximation, the propagat- 
ing equations for the external classical field amplitude 
E+ = E+ + E+ read 

id z n c + ^-V 2 T n c = -^|.g c | 2 a(° e >(22a) 

id t n p + icd z n p + ^-v 2 T n p = -n\ 9p \ 2 a^\22b) 

£l\l p 

where V 2 ^ is the transverse laplacian and 



9j 



= d, 



2e hV 



, J = c,p 



(23) 



The right hands of Eqs. (j22|) are the response functions 
of the atoms located in position r. 

For simplification, we deal with our problem in two- 
dimension, that is in the x-z plane. Such a system can 
be realized experimentally by putting the atomic ensem- 
ble in a planar waveguide. In the large detuning between 
the control light and the atomic transition, most atoms 
are in the ground state \g), which means the population 



difference in thermal equilibrium w e 



-1. As the de- 



tuning usually is much larger than the dephasing rate 72, 
we can neglect the imaginary part of the atomic response 
to the optical fields, then the propagating equations read 



td z n c + — d 2 n c = : (24a) 



l + 2/?|fl c 



idttip + icd z n p + ^-dln p = apflp 2 (24b) 

2kp (i + 2p\n c \ 2 



where 8 = has been set and 
N, ,, 



p = 



272 

7i (7 2 2 + A 2 ) 



(25a) 
(25b) 
(25c) 



The consideration for taking 5 = reflects the principle 
of this problem by the following reasons: 1) a resonance 
occurs at S = which can be obtained by the inspection 
of Eq. (|2"0|); 2) Though two-level atoms oscillate at fre- 
quency 6, it becomes important only when 6 is less than 
71, which results in the occurrence of a narrow spectral 
hole with width proportional to 71 burning in the homo- 
geneous absorption profile of the probe field. Hence the 
group velocity of the probe beam is somewhat reduced 
by CPO. 

As the atomic response to the probe beam depends 
on the control beam, we first deal with the propagating 
equation for the control beam. In the limit of low inten- 
sities of the control beam, the control wave experiences 
a refractive index correlated with the optical intensity 



id z Q e + ^d 2 x n c = oo(i- 2(3 |r> c | 2 ) n c , 



(26) 



i.e. the control wave induces a purely third-order nonlin- 
ear optical response. By defining the momentum opera- 
tors p M = — id^ and the effective mass m p = k, Eq. (|26p is 
a nonlinear Schrodinger-like equation, which has a soliton 
solution 



0, 



2 



sech ( — 
Lie 



(27) 



with Q = u c fl. Such solitary wave can propagate for long 
distances with an invariant transverse profile. Here 



L r = 



V2/kc 

Otcfi 



(28) 



is the transverse size of the control beam. 

We now consider the propagation of an initial Gaussian 
probe wavepacket 



Vt p (0,x,z) = 



1 



: exp 



(x - a) + z 2 



(29) 



in this two-level medium, where b (< L c ) is the width 
of the probe field and a is the initial wave packet center 
of the probe field along x-direction. The magnitude \a\ 
denotes the distance from the center of the control beam. 
In order to investigate the evolution of this initial state, 

we expand ^1 + 2(3 |fi c | 2 ^ at the position a and retain 
the linear term proportional to x — a. Then the paraxial 
equation for the probe beam becomes 



idt$l p + icd z Q p 



where parameters 



Vo 



—d 2 x n p = [ m + m {x~a)]n p (30) 



^sech 2 (— 
2 U c 



1 -2 



m 



a p Q 2 !3. 



sech 2 (f-) tanh [f- 
1 + ^sech 2 (f-) 



(31a) 



(31b) 
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Also by defining the momentum operators p^ = — id^ 
and the effective mass m p — k p /c, the dynamic of the 
probe wave fi p is governed by the Schrodinger-like equa- 
tion idt^lp = Hflp with the effective Hamiltonian 



H = cp z + - — pi +r) +r) 1 (x-a). 
2m„ 



(32) 



Thus by making use of the Wei-Norman algebraic 
method [22j, |25j (see the appendix), it can be found that, 
after passing through the Rb gas cell, the initial center 
(x, z) — (a, 0) of the probe field is shifted to 



x = a - 
z = L 



2kc 



(33a) 
(33b) 



By tracking the center motion of the probe beam, a mi- 
rage effect might occur to the probe light, which is de- 
termined by the sign of the detuning A and the initial 
center a of the probe beam. When the center is collinear 
to that of the control field a = xq = 0, the trajectory of 
the signal light is a straight line. We assign the positive 
sign for a as the probe beam shifted to the right with 
respect to the center of control light, and denote a < 
as the signal beam shifted to the left. When the probe 
beam is shifted to the right, i.e. a > 0, for red detun- 
ing A < 0, the probe wave feels a "attractive potential" 
within the atomic medium due to rji > 0, thus the tra- 
jectory bends into the left side of the e z -axis; for blue 
detuning A > 0, the signal light undergo an "repulsion 
potential" due to r/i < 0, therefore the ray bends to the 
right. When the probe beam is shifted to the left, i.e. 
a < 0, as the coefficient of the linear potential is larger 
than zero, that is A > 0, the probe beam experiences a 
"attractive force" within the two-level atomic medium, 
and its center is shifted to the left. As rji is smaller than 
zero, i.e. A < 0, the probe beam suffers a "repulsion 
force" during passing through the atomic medium, hence 
the whole wave packet goes to the right. The correspond- 
ing schematic diagram is given in FigfJJ where dark thick 
line is the transverse spatial profile of the control light. 
The dash lines give the deflection at A < 0, the dotted 
lines describe the light trajectory at A > and the black 
solid depict the light ray at a = 0. 



V. CONCLUSIONS 

In conclusion, we have theoretically predicted a phe- 
nomenon of the enhanced light deflection by an atomic 
ensemble through coherent population oscillation mech- 
anism, which is realized by a two-level atomic ensem- 
ble interacting with a control field and a much weaker 
probe field. After calculating the atomic response to the 
external fields by the perturbation approach, we obtain 
a system of equations that describe the spatial motion 
of the control field and the probe field in the two-level 
atomic medium. The propagation of the control field is 
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FIG. 2: Schematic illustration about the ray deflection of 
the probe light in the presence of spatial-distributed coupling 
light. The dark thick line is the transverse spatial profile of 
the control light. 



governed by a nonlinear Schrodinger equation in the limit 
of a low intensity. It shows that a solitary wave can be 
excited, which has an invariant transverse profile for a 
long-distance propagation. Due to a transparency win- 
dow with a width of the order of 71 and the achievement 
of the substantially reduced group velocity for the probe 
field by CPO, we found that the deflection of the light 
ray can be controlled by two controllable external param- 
eters: the initial center of the probe beam with respect to 
the control light, and the detuning between the control 
field and the atomic transition, which is similar to that 
in EIT medium investigated by Ref . [13, HH . Our anal- 
ysis maybe provide a new technique for transverse light 
guiding and it is much appreciated that our prediction 
can be verified in the further experiment. 
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APPENDIX A: FACTORIZATION OF UNITARY 
OPERATOR 

The shifted distance along the transverse direction in 
Eq. (|33p is obtained by applying the unitary operator 
U (t) to the initial wavefunction of the probe field, where 
the unitary operator U (t) is generated by Hamiltonian H 
in Eq. (|32[) ■ As operator p z commutes with other terms 
in Eq. (|32p , the unitary operator U (t) can be firstly 
factorized as 



where 

is generated by 

Hrr 



U{t)=e- m ° tp *U'{t) 1 (34) 



U' (t) = exp (-iH x t) (35) 



-Pi +Vo+Vi(x- a). (36) 



2m 
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It means that the unitary operator U' (t) only contains 
operator P x and x, which generate the Lie algebra with 
the basis {x, P x , P x , 1}. Thus, operator can be fac- 
torized as the form 



This leads to a systems of coupled differential equations 



U'{t) 



(37) 



and gi = gi(t) are unknown functions of time t to be 
determined. Here x' — x — a. 

Mathematically, the above factorization Ansatz is 
based on the Wei-Norman algebraic theorem [25|: if the 
Hamiltonian of a quantum system 



H 



K 

E 



(38) 



is a linear combination of the operators Xj that can gen- 
erate a iV-dimensional Lie algebra with the basis: 



.dgi 

.dg 3 
l ~dt 



1 

2m 



dgi dg 3 
dg4 .dg 3 

-dT ~ l ~dt 92 



= o, 



-Vo. 



(42) 
(43) 
(44) 

(45) 



{Xi,X 2 , Xk, Xk_i, ....Xn}, 



(39) 



then the evolution operator governed by H can be fac- 
torized as a product of the single parameter subgroups, 
that is , 



N 



U = l[e 



(40) 



3-1 



The solution to these equations reads 



where the coefficients £j (t) can be determined by the "ex- 
ternal field parameters" Cj (t) through a system of non- 
linear equations. 

Now, we differentiate (|3"T1) with respect to t and mul- 
tiply the resulting expression on the right hand side by 
the inverse of (|3"T|) . obtaining 



1 



2m, 



-Px + Vo + Vix 



(41a) 



dt * + x I dt m dt 



■ dg 3 
lX ~dt 



dg^ 
dt 



dt 



9i 
93 

92 
.94 



t 

~ l 2m' 
-irjit, 

2m ' 



Vot 



3 2m 



(46a) 
(46b) 
(46c) 

(46d) 



By applying the unitary operator U (t) to the Gaussian 
probe wave packet in Eq. (f29|) . a straightforward calcu- 
lation shows that the center of the probe wave packet is 
shifted to the position given in Eqs. (f33|) . 
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